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SOURCE - TEE RAYLEIGH CRITERION 
By Boa^-Teh Chu 

SUMMARY 


The stability of systems containing a heat source is examined from 
the energy point of view. Rayleigh's criterion is derived. In the case 
of a flame, it is found that Rayleigh's criterion must be modified slightly 
if the specific-heat ratios y of the burned and unburned gases are 
different. 


INTRODUCTION 


It has long been known that oscillatory phenomena are of common 
occurrence in systems containing an energy source, whether it is distrib- 
uted in space or concentrated in a limited region (e.g., in a plane or 
a surface) . Some typical examples are the screaming and chugging of 
afterburners and combustion chambers, as well as the more classical 
examples, such as the Rljke tone, the singing flame, and the oscillatory 
flame propagation in a tube . Generally speaking, the phenomena possess 
two characteristics. First, an oscillation is built up without any appre- 
ciable external excitation. Secondly, after the oscillation has increased 
to a certain amplitude, it is maintained in this state (unless the system 
breaks down before this state is attained) . 

A dynamical Bystem will start to oscillate with increasing amplitude 
only if energy is fed into the system in such a way that there is a net 
increase of the total mechanical energy of the system after each cycle 
of osc ill ation. The vibration is finally maintained at a given level 
when the mechanical energy fed into the system per cycle of oscillation 
is just equal to the sum of that dissipated by viscosity and that radiated 
away from the system per cycle of oscillation. When the amount of energy 
released per second is controlled by an external agency and is independ- 
ent of the fluctuation inside the system, .the oscillation will build up 
when the energy is released at certain characteristic frequencies. Such 
a phenomenon is usually described as resonance. On the other hand, if 
the system itself contains an energy source with the property that the 
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amount of energy release depends upon the fluctuation inside the system, 
an accidental small disturbance inside the system may interact with the 
energy source in such a manner that mechanical energy will be fed into 
the disturbance per cycle of oscillation, building up its amplitude to 
an appreciable extent. Depending on the nature of the response of the 
energy source to the disturbance, the interaction of the two may lead to 
a buildup or damping of the oscillation in the system. When the oscil- 
lation does build up, the phenomenon is visually described as instability. 

It Is the purpose of this analysis to study the conditions under which 
systems containing heat sources may be expected to exhibit instability. 

To illustrate more precisely what has been described and also as an 
Introduction to the problems proposed for study, consider a special 
example. Consider a tube of length L with open ends. Let the origin 
of the coordinate system be chosen at one end of the tube with the x-axis 
parallel to the wall of the tube. Let Q(x,t) be the rate of heat release 
to the medium per unit volume per unit time at the point x and at the 
Instant t. Let Cp be the specific heat at constant pressure of the 

gas in the tube and T 0 be the temperature of the gas when Q = 0. If 
Q /c p T 0 and its time derivative are "small," the pressure fluctuation In 
the tube will also be small and is governed by the equation 


_1_ dfp _ dfp _ _5_ / Q \ 
c2 3t 2 5x 2 at \ c p t °/ 


where P = Sp is the pressure fluctuation (above the mean pressure p Q ) , 

y is the specific- he at ratio, and c is the velocity of sound (see, 
e.g., ref. 1). At the two ends of the tube, there exists the condition 


P = 0 (x = 0, x = L) (lb) 


for all- values of t. If Q(x,t) is a given function of x and t, 
the norihomogeneous term at the right of the equation is known and repre- 
sents the "forcing function" or "source distribution." The problem is 
then one of forced oscillation. If the Initial conditions are specified, 
for example. 


(P)t=0 - f lW 



= ^ 2 ( x ) 




) 


(lc) 
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a solution can be constructed readily. It can be easily shown that when 
Q(x,t) is periodic in t with a period conmeasurable with the natural 
period of the tube, that is, 2l/c, the amplitude of a mode of oscilla- 
tion will increase linearly with time. This is the resonance phenomenon. 
Sow, if the rate of heat release Q is not a prescribed function of x 
and t but is given In terms of the fluctuations in the system., for 
example, Q is directly proportional to the pressure fluctuation P, 
the differential equation governing the pressure field becomes a homo- 
geneous equation. Thus, if 


Q, 

G p-o 


= kP 


(2a) 


then equation (la) becomes 


_1_ dfp cjfp = k SP 
c 2 St 2 ' Sx 2 " at 


(2b) 


The solution of this equation satisfying boundary condition (lb) and 
initial conditions (ic) is 


P = 


c 2 !^ 
e ^ 

n=l 



(An cos o^t + Bft 


sin o^tjsin 


ttrcx 

L 


(2c) 


where 


Oh » 



k 2 ^ 

~T~ 


(2d) 


An - 



(ae) 
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Bn = 


L<% J 0 


/ 
d n 


f 2 (x) - f]_(x) 


Bin — dx 


(2f) 


Prom equation (2c) it is seen that, if k > 0, the initial disturbances 
will grow exponentially with the time t and such a system is unstable. 
If k < 0, the solution clearly shows that the system is stable to small 

disturbances . 

The same conclusion can be arrived at more simply by using the for- 
mal method of stability analysis . As usual a disturbance of the form 

P = e^P’k-ijr (x) is assumed. To satisfy the differential equation it is 
necessary to have 

ijr( x ) = A cos A(p)x-+ B sin A(p)x (2g) 


where A and B are constants and 


A(p) 



(2h) 


To satisfy the boundary condition (lb) , It is necessary that A = 0 

and A(p) = n being a positive integer. Solving for p from the 

L 

second condition gives 



(2J) 


Consequently, if k > 0, both roots have the property that Im p n <0; 
if k = 0, Im p n = 07 and if k <0, both ro ots are such tha t Im p n > 0 
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where Im Indicates imaginary roots. Since P = eip' i> i|r(x), the system 
is unstable for k > 0, but stable for k fs 0. In this formal analysis, 
usually no time is taken to satisfy the initial conditions (lc) in the 
belief that these initial conditions can always be satisfied by the super- 
position of the assumed form of disturbances. In the present analysis, 
this is indeed the case since the eigenfunctions t|t(x) corresponding to 
n = 1, 2, . . . form a complete orthonormal set in the interval 
0 x L. 

In the above example, when the system is unstable (i.e., k > 0) , 
the amplitude of the fluctuation grows exponential! y with time without 
ultimately settling down into a state of permanent oscillation of con- 
stant amplitude (i.e., into a limit cycle). If, however, the rate of 
heat release Q is related to the pressure fluctuation by 


. i PS) (3a) 

where e and P Q are two positive quantities, instead of by equa- 
tion (2a), the differential equation governing the pressure fluctuation 
becomes 


i tZ . & = Jp 0 2 . p2)|E 

OO O VO /£ t 


e 2 St 2 


bx 


(yo) 


Any small disturbances will at first be amplified exponentially and will 
ultimately settle down into a limit cycle since e(P Q ^ - P^) is positive 
for jpj < P Q and becomes negative for |pj > P D . This is a case of 
selfsustained vibration of large amplitude. This particular system clearly 
corresponds to the limit of a very large number of identical Van der Pol 
oscillators linearly coupled together to form a continuum. It may be of 
some basic interest to study more closely a nonlinear system like 
equation (jb). 

In the above examples, heat is released at all points in the tube. 

If heat Is released in a portion of the tube or concentrated in a narrow 
region, similar analysis applies. However, the calculation is, as a rule, 
longer a.nd the simplest method of attack is probably the formal stability 
analysis. However, with each slight change in the geometrical arrangement, 
for example, the extent and connectivity of the region in which heat is 
released, a new calculation must be performed. In this sense, the above 
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analytic approach is not very satisfactory although it does provide 
quantitative information as to the rate of amplification of the disturb- 
ance. For many purposes, however, precise quantitative information is 
not so important as having a qualitative idea as to whether a system may 
or may not be stable and why it is so. The latter is indeed the more 
important if the fundamental physical principles involved are to be under- 
stood. It is toward this goal that the study in the remainder of this 
paper is directed. 

The present investigation was conducted at the Department of 
Aeronautics of The Johns Hopkins University under the sponsorship and 
with the financial assistance of the National Advisory Committee for 
Aeronautics . 


SYMBOLS 


A 

A 

A f 

An 

B 

Bn 

C P 

c 

E 

fl(x),f2(x) 

h 

k 

L 

n 

P 


cross-sectional area of tube 

constant 

flame area 

defined by equation ( 2 e) 

constant - 

defined by equation ( 2 f) 
specific heat at constant pressure 
velocity of sound 

total energy. Kinetic energy + Potential energy 

arbitrary functions 

location of heater 

proportionality constant 

length of tube . 

positive integer 

pressure fluctuation above mean level 
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Po 


P 

Q,Q(x,t) 

Q 

R 

s a 

St 

T 

To 

t 


u 

X 

7 

A 

S( ) 
€ 

x(p) 

P 

Hr(x) 

C13 


o>n 


positive quantity 
pressure 

heat release to medium per unit volume per unit time 
heating value of mixture 
gas constant 
apparent flame speed 

flame speed 
temperature 

temperature when Q = 0 
time 

velocity 

space coordinate 

parameter 

specific-heat ratio 

increment of total energy after one cycle of oscillation 

increment of ( ) 

position quantity 

defined by equation (2h) 

density 

eigenmode 

rate of heat release per unit cross-sectional, area of tube 
value of co when there is no disturbance 
defined by equation (2d) 
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Subscripts: 

1 conditions upstream of heat source 

2 conditions downstream of heat source 

h conditions at heat source 

i conditions at inlet end of tube 


RAYLEIGH ' S CRITERION 


In explaining the process of maintenance of vibration by heat, 

Rayleigh (ref. 2, pp. 226 and 227) states that if heat is periodically 
communicated to and abstracted from a mass of air vibrating in a cylinder 
bounded by a piston, for example, the effect produced will depend upon 
the phase of vibration at which the transfer of heat takes place. If 
heat is given to the air at the moment of greatest condensation or Is 
taken from it at the moment of greatest rarefaction the vibration is 
encouraged. On the other hand, if heat is given at the moment of greatest 
rarefaction or abstracted at the greatest . condensation the vibration Is 
discouraged; however, there is no effect in encouraging or discouraging 
the vibration if the air concerned is at a loop, that is, a place where 
the density does not vary, or if the communication of heat is the same 
at any stage of rarefaction as at the corresponding stage of condensation. 
Rayleigh went on in applying this criterion to explain the singing flame 
and Rijke's tone. He was able to give a very satisfactory account of the 
observed effect of tube length on the excitation of singing flame and the 
effect of the position of the heater in the Rijke's tube in the production 
of Rijke tone. (See pp. 229 to 230 and 232 to 233 of ref. 2.) 

In most later investigations on the maintenance of vibration by heat, 
Rayleigh's criterion has been quoted (e.g., refs. 3 to 7)* Rayleigh him- 
self did not show how he arrived at such an apparently rather general 
statement. Noting this deficiency, Putnam and Dennis proposed a proof 
which many find hard to believe (see appendix of ref. 3) • Their proof 
is not wrong, but it contains so many plausible hypotheses of a mathe- 
matical nature that one Is not sure if the physical conclusion it purports 
to prove is not more plausible than the mathematical hypotheses used. 
However, Putnam and Dermis did put Rayleigh's criterion In a very precise 
form: Mathematically, and neglecting damping forces, Rayleigh's criterion 

requires that the time integral, over a cycle, of the instantaneous prod- 
uct of the rate of heat release and of the oscillating component of the 
pressure be greater than zero (ref. 5)-^ 

^The criterion in this form was also independently suggested by Bailey 
(see p. III . 9 of ref. 6). 
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An explanation of Rayleigh. 's criterion on a physical "basis will 
make it quite cl ear why the criterion is in fact a very plausible one. 2 
For convenience of explanation, consider a tube of finite length. Divide 
the tube into many fictitious compartments (say of equal size) in one of 
which heat is being added periodically from an external source. This is 
illustrated in figure 1 where the compartments are separated by planes 
(shown as dashed lines in the figure), while the shaded compartment bounded 
by the planes marked A and A' is the one in which heat is being added. 

In the absence of viscosity and heat conductivity, there is no loss of 
generality by replacing the fictitious planes A and A’ by two solid plane 
surfaces (or "pistons”) provided that these walls always move in synchro- 
nism with the motion of the fictitious planes A and A 1 . The other com- 
partments serve, as is well known, as a mechanical spring -mass system 
(or flywheels in the engine analogy) where mechanical energy produced by 
the "engine" compartment can be stored. The work done by the engine can 
be easily calculated from the FV -diagram. Depending on when heat is 
added and subtracted during each cycle of vibration the amount of mechan- 
ical work produced is greater than, equal to, or less than zero; this 
determines whether mechanical energy is given to or extracted from the 
spring -mass system (i.e., the flywheels). For the case when the amount 
of heat energy released by the heat source depends on the magnitude of 
fluctuation in the system and is zero when there is absolutely no fluctua- 
tion in the system, the initial disturbances act as "an "engine starter." 

It is at once clear from the FV-diagram that if heat is added when the 
pressure in the engine is the highest and taken away when the pressure 
in the engine is the lowest, the maximum amount of mechanical work is 
obtained. This is in fact identical to Rayleigh's statement that if 
heat is given to the air at the moment of greatest condensation or Is 
taken from it at the moment of greatest rarefaction the vibration Is 
encouraged. Furthermore, it is also clear from the engine analogy that 
the real criterion of amplification of a disturbance is that the net 
mechanical work done by the engine per cycle must be greater than the 
loss through viscous dissipation and, hence, must at least be greater 
than zero. This is also in accordance with Rayleigh's remark that there 
is no effect (of encouraging or discouraging the vibration) If the air 
concerned Is at a loop, that Is, a place where the density does not vary, 
or if the communication of heat Is the same at any stage of rarefaction 
as at the corresponding stage of condensation. 

It is possible to put the above criterion in a mathematical form. 

The result of this calculation leads to the formulation by Putnam and 
Dennis of Rayleigh's criterion. Since this calculation will appear as 
a particular case of the analysis of a more useful system discussed in 
the next section, the details will not be given here. 


p 

The author Is indebted to Prof. Leslie S. G. Kovasznay for this 
explanation. 
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In the example Just presented, heat must he guided as veil as sub- 
tracted from the "engine compartment" if the state of the gas inside the 
compartment is to he brought back to its initial state after each cycle. 
The heat subtracted is, of course, less than the heat added if the engine 
produces a net positive, amount of mechanical vprk every cycle. The heat 
subtracted corresponds to. the heat rejected in an actual engine. To be 
sure, for most practical cases, the gas in the engine compartment instead 
of returning to its initial state at the end of a cycle ends up at a tem- 
perature slightly higher than its original temperature . This cuts down 
the mechanical output of the engine in the next cycle and sometimes also 
reduces the amount of heat given to the engine (e.g., vhen the rate of 
heat transfer depends on the temperature of the medium as in the case, of 
conductive heat transfer) . These are precisely some of the factors vhich 
vill eventually limit the amplitude of the oscillation to finite size and 
account for the self sustained vibration (i.e., limit cycle) vhich is often 
observed in experiments on such systems . 

When heat is not added to the gas in the limited region as Bhovn in 
figure 1 but is distributed throughout the tube, the same argument applies 
except that it is necessary to think in terms of a multicylinder engine 
rather than a single -cylinder engine doing vork on a spring -mass system. 

In fact, vhen the amplitude of the oscillations is small this case can 
alvays be analyzed by a superposition of the previous case. 

The only drawback of the engine analogy is that it cannot be extended 
simply to the case vhen there is a constant current of air through the 
tube and vhen the heat added to the medium fluctuates with the disturbance 
in the system above some nonzero mean value. It is clear that the heat 
added to the system can be decomposed into tvo parts, a steady component 
(i.e., a direct -current component) vhich maintains a given tenperature 
distribution inside the system and a varying component (i.e., an 
alternating-current component) which feeds energy into the disturbance. 

The heat energy that has been given to the disturbance appears in two 
forms: The pressure waves which carry the mechanical energy into the 

rest of the system, and the entropy spots vhich are carried away by the 
draft and represent the amount of energy going into the heating of the 
gas above the mean temperature distribution. The latter is therefore 
unavailable for doing mechanical work. In such a case, the entropy spots 
carried away by the draft vill be the counterpart of the heat energy 
rejected by the engine. It is possible to study the energy transforma- 
tion and, in particular, the total amount of the mechanical energy stored 
in the system per cycle of vibration from a purely analytical standpoint 
(i.e., without appealing to a physical model). This will be done in the 
following sections. The physical model, however, is useful in inter- 
preting the analytical formulation and conclusions. 
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ANALYTICAL FORMULATION 


It is the purpose of the present analysis to establish the condition 
under which a small disturbance is amplified in the course of time by its 
interaction with a heat source. Though the analysis is restricted to 
small disturbances for analytical, reasons, the derived condition is never- 
theless useful in practice since small disturbances always exist in any 
physical setup and also because in most cases a large disturbance does 
not disappear without first becoming small. 

For definiteness, consider a gaseous medium in a tube of length L. 
The origin of the coordinate axis is chosen at one end of the tube and 
the x-axis, parallel to the wall of the tube. The two ends of the tube 
will then be given by x = 0 and x = L. The heat source is assumed to 
be concentrated in a single plane at x = h .3 It is further assumed that 
there is a constant current of air through the tube flowing into the tube 
from the end x = 0 with a speed much smaller than the local speed of 
sound. Finally, it is assumed that in the absence of any disturbance 
the heat source is releasing energy at a rate of a> 0 units of heat per 

second per unit of cross-sectional area. As a result of the heat addi- 
tion, the state of flow ahead of and behind the heater will not be the 
same. Let pij_, p^, Tj_, and u-j_ denote, respectively, the pressure, 

density, temperature, and velocity of the undisturbed medium upstream 
of the heat source (i.e., for 0 < x< h) , while Pg, Pg, Tg, and Ug 

denote, respectively, those downstream of the heat source (i.e., for 
h < x < L; see fig. 2 ) . These quantities are related by the continuity, 
momentum, and energy equations across the heater which are 

p l u l = ^^2. (to.) 

p l + p l u l 2 = p 2 + P 2 U 2 2 C 1 *) 


“O - + i “2 s ) - Pi^ppA + \ Pi 2 ) 


In case the heat source is a flame front which propagates with a flame 
speed relative to and against the oncoming stream, the heat source 


5 This assumption by no means Limits the generality of the results 
obtained below since, for small disturbances, a continuously distributed 
heat source can always be regarded as a superposition of a large number 
of concentrated heat sources. 
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will be seen to have an apparent velocity S a (taken as positive if the 
heat source moves forward against the oncoming stream) given by 

s a = St " U 1 (5a) 

In such case, equations (4) must be replaced by the following equations: 

p l( u l + s a) = P 2 .(^ 2 . + s a) (5b) 

Pl + p-L^ + S a ) 2 = p 2 + P 2 (u 2 + S a ) 2 (5c) 

03 o = P 2 (u 2 + S a ) c P 2 T 2 + |(u 2 + S a ) 2 - P!(u x + S a ) jcppTp + |(u x + S a 

(5^) 

Evidently, the system given by equations (4) can be considered as a par- 
ticular case of the system given by equations ( 5 ); for, if = u x in 

equation (5 a ) it is found that S a = 0 and equations ( 5 ) reduce to equa- 
tions (4) . Consequently, only system ( 5 ) needs to be considered. 

If p^, p-^, Tp, u^, o> 0 , and are assumed to be given, then 

equations (5a) to (5d), together with the gas law 

p 2 = p 2 R 2 T 2 (5e) 

will enable calculation of the five unknowns p 2 , p 2 , T 2 , Ug, and S a . 

In particular, it can be easily verified 
the same order of magnitude as ^u x + S a 

velocity of the direct current through the tube and are small com- 

pared with the local sound speed (so that terms which are of the order 
of the square of the Mach numbers can be dropped), equations ( 5 ) can be 
reduced to 


■t (u 2 + Sa) j \J 7 2 R 2 T 2 1s ^ 
7 ^R^T^. Consequently, if the 
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Pl( u l + s a) = P 2 ( u 2 + S a) 


(6a) 


Pi = P 2 


(6b) 


0)0 




?b 

7i - 1 


P l( U l + 



(6c) 


where 


s t = u i + s a 


(6d) 


P 2 — 


(6e) 


Now, suppose that there are some disturbances inside the tube. Let 
the pressure, density, temperature, and velocity of the medium upstream 
of the heat source at the point x and instant t be denoted by 
p-|_ + Sp-j_> Pp + SPi> Ti +■ and u-]_ + 5u^_ while those behind the 


heat source are denoted by P 2 + SP 2 , 
If the disturbances are weak (l.e., if 


P2 + Sp 2 , ^2 + ^2 + 5 u 2* 

Spi Spi STi Su^ 5p2 

~P? Pi T L ’ ~ 2 3 


8P2 Sip &Uo \ 

, , and « 1, where c denotes the velocity of sound), the 

P 2 I 2 c 2 / 


energy in the disturbances can be decomposed into two parts: The kinetic 

energy is given by 


2 P 1 ( 5u l) 2A ^ + 2 P 2 f h ( 5u 2 ) 2a ^ (7a) 


(A being the cross-sectional area of the tube) and the potential energy 
is given by 
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(see, e.g., pp. 17 and 18 of ref. 2 ). The total energy in the disturbances 
is the sum of these two equations. If the total energy in the disturbances 
increases after each cycle of oscillation (described subsequently) , it is 
said that the disturbance is amplified by the heat source. Now, it is a 
well-known fact that the surface integral of the product of the pressure 
and the velocity component in the direction of the normal of a control 
surface Is related to the rate of change of energy inside the system. 

(See, e.g., ref. 8 .) This suggests the following calculation. 

Consider a control surface shown by dashed lines in figure 3* It 
consists of two compartments : The first encloses the medium ahead of 

the heat source, and the second encloses the medium behind the heater. 
Consider the first compartment. If the subscripts h and i are used 
to denote the condition of flow at the heat source and at the inlet end 
of the tube, respectively, the following identities are obtained: 


(pi + 8Px) h ( u i + ^i) h - (pi + & Pl)l( u l + Su i)i = f Q ^[(*>1 + & Pi)( u i + 6u i)] 

= f 0 [( Pi + 8P1 ) b ( Ui + 6ui ) + 

(^1 + Su;l) ~ (pi + 8Pi)l a* 

= Pi^5u^ h - 


JC[W h K> * W b H] 


dx 


that is, 


(&Pl)h( 5 u l)h - ( 5 Pl)i( 5 u l)i = j Q ^Pi |- (sup) + Sup |- (& P p) 


dx 


( 8 a) 
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When the I nl e t end Is open, (§Pi}i = 0; hence. 


(Sp 1 ) h (8u 1 ) h - ^ [tel (BUi) + 8U X |- (5pi) 


dx 


(8b) 


(This formula is still true if the inlet end is closed since (dui)^ = 0.) 

Likewise, if the compartment of the control surface behind' the heater 
is considered. 


( 8P2 )h( Su 2)tL = -f ( 6p a) f; ( 5u 2) + ( Su s) £ ( 6 ®a) (8o) 


It follows from the last two formulas, equations (8b) and (8c), that 


rhp 

( 5 ®2)h( Su 2)h " (pPl)h( 5u l)h = ” J B Pl ( 5u l) + Bu l ( B Pl) 


dx - 


i [ 6P2 b ( SU2 ) + 5us b ( SP2 ) 


dx 


( 84 ) 


Note that this equation is exact for tubes with open (or closed) ends; 
that is, no neglection has been made. Next, the simplifying assumption 
that the disturbances in the system, are small is introduced into equa- 
tion (8d). It will be shown that the right-hand side of the above 
expression gives precisely the rate of increase of total energy Inside 
the tube if the Mach number of the steady current through the tube is 
of the order of magnitude of 0.01 or less fe.g., if uj_ < 10 fps, and 

Ct = 1,000 fps, then ^ ^ 0.01). For such low Mach numbers, the hydro - 

C 1 / 

dynamic equations governing the disturbance can be written simply as 




= 0 


a_ 

3t 


(9a) 
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S_ 

St 



+ c 2 


S_ / 5p l \ 
Sx V^P-^/ 


= 0 


(9b) 


S_ 




(9c) 


8T]_ 8p^_ 8p]_ 

T 1 Pi Pi 


(94) 


for the region ahead of the heat source. Likewise, similar equations 
(with subscript 2 instead of 1 ) hold for the region behind the heat 
source. 

Consider the first integral on the right-hand side of, equation ( 8 d) 
The following equation is obtained by use of equations ( 9 ) : 


1 [ 8Pi b ( 5Ui ) + SUi b ( SPi )] " f 0 


s / 5p l 

5p i f- + 

St WqP-j 


7T W b M 


dx 


S_ 

St 



1 ( & pl ) 2 , 

2 Vi 



( 10 a) 


and gives the rate of increase of the disturbance energy (per unit cross 
sectional area of the tube) in the region ahead of the heat source. 
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Likewise , 



a_ 

dx 


( Bu a) + Su 2 b ( EP2 ) 


dx 


6_ 

St 



1 

2 



dx 


(10b) 


gives the rate of increase of the disturbance energy (per unit cross- 
sectional area) in the region behind the heat source. Substituting equa- 
tions (10a) and (10b) into equation (8d) leads to the following conclusion: 
The rate of increase in the total energy E in the disturbance is given 
by 


b (E) “ A [( 6S2 ) 11 (. 5ua )K - KM 


(11) 


Next, the right-hand side of the above equation must be evaluated in 
terms of the rate of heat release at the heat source. Since equations (6) 
are the mathematical representations of the conservation laws at the heat 
source, they must be satisfied at all instants provided that instantaneous 
values of the pressure, velocity, and other physical variables are used 
instead of pq, uq, and so forth. In particular, equations (6b) and (6c) 
must be satisfied; that is, 

Pj_ + 5p-j_ = P 2 + SP 2 (12a) 


ub + So = — — ^P2 + &P 2 ) ( u 2 + S a + 5^2 + 5S a ^ - 


7 l 


— — ^Pq + Spq^ (u-L + Sa + 5U]_ + 5S a ^ (l2b) 


In writing the second equation, allowance has been made for the change 
in the rate of heat release with the fluctuations at the heat source. 
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Since these relations are -valid only at the heater, the subscript h 
should have been used to indicate this fact. However, for the sake of 
simplicity, this subscript will not be attached to the various variables 
involved until the final result is obtained (i.e., eqs. ( 13 )). Making 
use of equations (6b) and (6c), the last two equations can be rewritten 
as 


Spi = 5p 2 

Sen 72 5u 2 + 6S a 7 ± Su l + ss a 

-... l . — + 

Pl c l 7 2 “ 1 C 1 7 X - 1 

72 S ^2 ( u 2 **• s a) + ( 5u 2 + ss a) 

72 “ 1 P 2 C 1 

7l Bpi (u x + S a ) + (suj + SS a ) 

7q " 1 Pi C 1 


(12c) 


(124) 


Equation (12b) has been put in a nondimensional form so that the order 
of magnitude of the various terms can be readily compared. If terms 

5pi &p 2 8 ut 8u 2 Un 

which are products of the small quantities , , , , — , 

Pl p 2 c-l c 2 c x 

u 2 k . - 

and — are dropped, equation (12d) becomes 


See _ 7 2 5^2 + 6S a 5u x + 8S & 

Pl c l = 7 2 “ 1 c-j_ “7i-l c-l 


(12b) 


^ote that this is the same assumption which was used in deriving 
equations (9)* 
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Introducing now the subscript h to indicate that these relations are 
valid only at the heat source, equations (12c) and (l2e) can be rewritten 
in the form 

(6p 2 ) h = (SPi)^ (13a) 


(-4 - 0*4 SSa >* 

(13b) 

Substituting equations (13) into equation (ll) gives the important 
relation 


~ W - A(SP4 



7 2 P X 





+ 6s 4 


(i 4 ) 


Now, if the energy in the disturbance after each cycle of oscillation 
is examined, ^ the disturbance is said to be encouraged by the heat source 
if the total energy in the disturbance at the end of the cycle is greater 
than that at the beginning of the cycle. If the increase of the dis- 
turbance energy after each cycle of oscillation Is denoted by A, the 
following equation results : 



5 since the motion is not strictly periodic (i.e., the motion does 
not repeat itself exactly after a certain interval), the term "cycle of 
osci ll ation" must be clarified a little. It can be conveniently defined 
as the interval between two successive times at which (bPi)^ = 0. 
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When the various forms of losses in the system are negligible, the con- 
dition that a disturbance is amplified by the heat source in the course 
of time is then 


A > 0 


(16) 


If A < 0, small disturbances in the system are damped by its interaction 
with the heat source. If A = 0, the heat source neither encourages nor 
discourages any disturbances. (it merely modifies their wave form.) 

In the next section, the above results are applied to a few special 
cases. - 


APPLICATIOHS 


The results of the foregoing analysis can be applied immediately 
to two special cases of practical Interest: 

(1) When the heat source is a plane heater and there is a current 

of air through a tube with Mach numbers of the. order of magnitude of 0.01 
or less. 

(2) When the heat source is a flame front. 


Plane Heater 

In the case of a plane heater the ratios of the specific heat at 
constant pressure to the specific heat at constant volume of the medium 
ahead of and behind the heater may be taken as the same (i.e., 

71 = 72 = 7 ) • Equation ( 15 ) then becomes 


A = -A. f (5 Pl ) h (B<s) dt (17) 

7 Pi J Cycle v ' 


Disturbances in the tube will be amplified in the course of time if 
A> 0. This is Rayleigh's criterion in the form first suggested by Putnam 
and Dennis (ref. 3) • It states that if the time integral of the product 
of the pressure fluctuation and the fluctuation in the rate of heat 
release over a cycle of oscillation is greater than zero, the disturbance 
is amplified; if the integral is less than zero, the disturbance is damped; 
if the integral is zero, the disturbance Is neither amplified nor damped 
by the heater. 
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Since neither cco nor u^_ enters into equation (17)> the above 
conclusion applies equally well to the case in which = 0 
and/or up = 0 , that is, the case where there is no current of air 
through the tube and/or there is no steady component of heating. 


Flame Front 

In the case of the flame front, 7 ^ is, in general, different 
from. 7 2 • Furthermore, from the definition of the flame speed and equa- 
tion ( 6 d), 


+ 8S a = 5S^ 


( 18 ) 


Hence, 


A = 


72 ~ 1 
72 


A 


L 


Cycle 


(_s Pl ) h (sc) dt 

Pi 





~ 1 ^2 


4) a 4- (sp ^ 6st 


dt 


(19) 


This formula for the increase of total energy per cycle of oscillation 
is strictly true for a plane flame front propagating at a speed of the 
order of magnitude of 10 feet per second or less. Physically, it seems 
probable that It may perhaps also be applied to a nonplanar flame front 
when the axial extent occupied by the curved flame front is Bmall com- 
pared with the length of the tube, provided that 5a> is interpreted as 
the change of total heat release per second per cross-sectional area of 
the tube (see ref. 9 ) • Since the rate at which heat is released by a 
flame front (plane or nonplanar) is equal to the product of the flame 
area, the flame speed, the density of the fresh gas, and the heating 
value of the mixture, any changes (induced by the fluctuation in the 
system) of any of these variables will cause- a change in the rate of 
heat release per unit time per unit cross-sectional area of the tube. 

In fact, since no = Af S-^Qp^A, the effect of a change in the flame 


area Af and the flame speed on the rate of heat release is, to the 

order of approximation which has been used in deriving equations (l2e) 
and. (9), given by 
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5m 

E>l c l 


£t Q, 5A f + ££t Q £f 
C 1 R 1 T 1 A C 1 R 1 T 1 A 


( 20 ) 


The effect of changes In the heating value on 6o> has been ignored here 
mainly because it is normally very small. Substituting equation (20) 
into equation ( 19 ) gives 



where usejtias been made of equations (6a), (6c), and (6d) and 
= AfS-tQp^ A to establish the relation 


yg - 1 q Af 1 , y± y 2 - 1 £1 _ 1 

72 R l^l A 7i ~ 1 72 p 2 


( 22 ) 


employed in simplifying the result of the substitution. Consequently, 
if the pressure fluctuation is positively correlated® with the fluctua^- 
tions in the flame area and flame speed, vibration is encouraged. If 
they are negatively correlated, vibration is discouraged. 

In the study of the stability of systems containing a heat source, 
the basic problem is then to determine the response of the heat source 
to the disturbances. Detailed analysis of the dynamic characteristics 
of the heat source is usually very difficult. The major contribution 
of Carrier in the study of Rijke tone 1 Is precisely in this direction 
(see ref. 10) . When a detailed analysis for a particular system proves 
to be too complicated (as Is usually the case when a flame front Is 
involved) , the dynamic response of the heat source can often be assumed 
if the study is accompanied by experimental investigations . The excel- 
lent work of Bailey (ref. 6) Is a typical example of this approach and 
another example can be found In reference 7- Still another example in 
this connection can be constructed by making use of the experimental data 
of Kaskan (ref. 3). In his study of oscillatory flame propagation in 


^Two variables are said to be positively correlated if the time 
integral of their product over a cycle of oscillation is greater than 


zero. 
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tubes, Kaskan proved beyond any doubt that in such a case the vibration 
is amplified mainly because of the proper phasing of the change of the 
flame area with the pressure fluctuation. Ee went on to propose a mechr- 
anism which explains this interaction between the flame front and the 
pressure waves. Kaskan' s picture of the interaction is essentially as 
follows: When a flame extends itself into the oscillating boundary layer 

near the wall, the motion of the central portion of the flame front is 
much greater than that inside the boundary layer. Consequently, there 
is a change in the flame area whenever there is an acoustical oscilla- 
tion in the tube. This picture of the interaction is supported further 
by the experimental evidence that if the oscillatory frequency is so 
high and, hence, the boundary layer is so thin that the quen ching distance 
of the flame near the wall is greater than the boundary-layer thickness 
(i.e., the flame does not extend itself into the boundary layer), then 
there is no change in flame area with any acoustical oscillation inside 
the tube and there is no evidence that the small disturbances are ampli- 
fied and that the flame vibrates. The excellent experimental study of 
Kaskan should be followed up with a theoretical study of this phenomenon. 
The main question then is to calculate the flame configuration (and, hence, 
the accompanying rate of heat release) in an oscillating flow inside a 
tube - much in the same spirit of Carrier in his study of Rijke tone. 

The difficulty here is that even an approximate calculation of the flame 
configuration proves to be too complicated. However, a very crude esti- 
mation of the rate of heat release at the flame front can be made on the 
basis of Kaskan' s results. Thus, by mn.Td ng the two assumptions (l) The 
flame configuration remains essentially similar during its vibration so 
that the flame area at any instant is proportional to the height of the 
curved flame and (2) the central portion of the flame moves back and 
forth as if there were no boundary-layer effects while the edge of the 
flame in the boundary layer only moves forward with a constant velocity 
equal to the mean speed of propagation of the flame, it is easy to show 
that the change in flame area is directly proportional to the fluctuation 
in the displacement of the fluid particle at the flame front. This can 
serve as a starting point of a theoretical calculation of the oscillatory 
flame propagation in tubes along a line similar to that Bailey has fol- 
lowed in reference 6. 


The Johns Hopkins Hal vers it y, 

Baltimore, Md. , May 31, 1955- 
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Figure 1.- Sketch showing engine analogy. Dashed lines indicated planes 
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Figure 2.- Heat addition in a flowing stream. 
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Figure 3*- Control surface shown as dashed line (two compartments). 
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